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Abstract. The author’s previous theory of the artificial satellite is derived in terms of the dis- 
turbed eccentric anomaly. The present development, in terms of the orbital true longitude, is a 
substantial improvement over the earlier work in that it leads to the faster convergence for large 
eccentricities and to a smaller number of terms in the series representing the perturbations. More- 
over, each approximation of the radius vector and of the parameters determining the position of 
the orbit plane is obtained not in the form of a truncated infinite series but in the form of trigo- 
nometric polynomials in two arguments. These arguments are the mean true anomaly and the 
mean argument of the latitude. The present theory, like the previous one, permits the computation 
of perturbations of any desired order. Thus, any future information about earth’s gravitational 


field can easily be included. 

Introduction . In a numerical theory of arti- 
ficial satellites developed by the author [Musen, 
1959] the eccentric anomaly is used as the 
independent variable. Thus the major limitations 
of the earlier theory are its slow convergence and 
its unsuitability for orbits of large eccentricity. 
Other features of the method given in the 
previous article are the introduction of the 
disturbing function with the separated elliptic 
and nonelliptic anomalies, the appearance of 
the divisor Vl — e 0 * hi the equation determining 
the motion of the perigee, and a comparatively 
slow convergence of the series for the reciprocal 
of the radius vector. The divisor Vl — e 0 * can 
become troublesome in the case of elongated 
orbits. Moreover, the series for the reciprocal 
of the radius vector converges slowly, even for 
moderate values of the eccentricity. This slow 
convergence is reflected in the slow convergence 
of all following intermediate results. The present 
theory removes these difficulties to a substantial 
degree by basing the development on the use 
of the orbital true longitude v as reckoned in 
the osculating orbit plane from the departure 
point. With this fundamental modification the 
troublesome divisor Vl — e 0 * does not appear, 
and all consecutive approximations are obtained 
in the form of trigonometric polynomials and 
not in the form of truncated infinite series. In 
addition, there is no necessity to introduce the 
modified disturbing function with the two 
types of eccentric anomalies. 


The present method uses the idea of Hansen 
[1838] of the separation of all the perturbations 
into the perturbations in the orbit plane and 
the perturbations of the orbit plane, but the 
basic W function is slightly different from the 
classical form. The two basic arguments are the 
mean true anomaly cv — 7 r 0 and the mean argu- 
ment of the latitude gv — 0 O . 

All notations are the same as in the previous 
article: 

i = the inclination of the orbital plane to 
the earth’s equator. 

6 = the longitude of the ascending node. 
a = the angular distance of the node from 
the departure point. 

a, e,n— o _8/2 , the osculating semimajor axis, 
the eccentricity, and the mean motion, 
respectively. 

/ = the osculating true anomaly. 

Disturbing function. The disturbing function 
is taken in the same form as before : 

n = k 2 u( 1 - 3<P 2 ) + &3« 4 (3iA - 5f a ) 


+ Jb 4 ti a (3 — 30^ 2 + 35 V) 

a) 

where 


u — l/r 

(2) 

and 


\p = sin t sin ( v — <r) 

(3) 
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a = (l g)v + 6 0 — N — K (4) 

0 = (1 - h f )v + 6 0 - N + K (5) 

Xi = sin hi cos N, X 3 = cos \i sin K ^ 

X 2 = sin sin N , X 4 = cos ^ i cos K 

We deduce 

\p = 2(XiX 4 — X 2 X 3 ) sin — 0 o ) 

+ 2(X 2 X 4 + XjX 3 ) cos ( gv 0o) (7) 

and 

ip 2 = 2(X 3 2 + X 4 ‘')(X 1 2 + X 2 2 ) 

- 2[(X x 2 - X 2 2 )(X 4 2 - X 3 2 ) 

— 4X!X 2 X 3 X 4 ] cos (2 gv — 2 0 O ) 

“h 4 (XjX 4 X 2 X 3 )(X 2 X 4 

+ XjXa) sin (2 gv — 2 0 O ) . (8) 

= 6X! 3 X 4 3 sin (gv — 0 O ) 

— 2Xi 3 X 4 3 sin (3#t> — 30 o ) 

T- 6X, 2 X 4 2 (X 2 X 4 X]X 3 ) cos {gv 0 O ) 

— 6X t 2 X 4 2 (X 2 X 4 + XA) cos (3 gv - 3 0 o ) (9) 

= 6Xj 4 X 4 4 - 8X t 4 X 4 4 cos (2 gv - 20 o ) 

+ 2X x 4 X 4 4 cos (4 gv — 40 o ) 

+ 16Xi 3 X 4 3 (X 2 X 4 + XiXs) sin (2 gv — 20 o ) 

— 8X! 3 X 4 3 (X 2 X 4 + XiX 3 ) sin (4 gv — 40 o ) (10) 

The development of the perturbations can 
be done either in a completely numerical way, 
using the method of successive approximations, 
or in an analytical or semianalytical way. In 
the*method of successive approximations the 
following system of formulas can be used : 

~ = 3k 2 u 2 {l — 3 \p 2 ) + 4k 3 u{3\p — 5\p 3 ) 
du 

+ bk 4 u ( 3 - 30^ 2 + 35i£ 4 ) (11) 

= —6 k 2 u\p + 3k 3 u 2 {l — 5\p 2 ) 


J dtt 
u 2 d\ p 


In equation 14 the argument is taken in the 
form gv — 0 O = v — (1 — g)v — 0 O , and the 
differentiation is performed only with respect to v. 


dfi dti d\p \ 


1 d± 

2 dX x 

ldxP_ 

2 ax 2 


1 0\P 


1 

2 ax 4 


= +x 4 sin (gv — 0 O ) 

+ X 3 cos {gv — 0 O ) 
= — X 3 sin {gv — 0 O ) 

+ X 4 cos {gv — 0„) 
= “X 2 sin {gv — 0 o ) 

+ X x cos {gv — 0 O ) 

= +X x sin {gv — 6 0 ) 

+ X 2 cos {gv — 6 0 ) 


(15) 

(16) 

(160 

(16") 


(13) 


+ 2Qhu{-3Js + 7<A 3 ) (12) 

J. dft 1_ dfi chp 

u dv u d\p dv 

\hr = (*>*« — ***») cos ( p* ~ e ») 

2 dv 

— (X 2 X 4 + X,X 3 ) sin {gv — 6 0 ) (14) 


(16"') 

If we wish to have an analytical or partly 
analytical development, we obtain 

= k 2 [l — 6(X x 2 + X 2 2 )(X 3 2 + X 4 2 )]u 3 
* + 6/b 2 [(X , 2 - X 2 2 )(X 4 2 - X 3 2 ) 

- 4XjX 2 X 3 X 4 ]?/ 3 cos (2 gv — 20 o ) 

12A; 2 (XiX 4 X 2 X 3 )(X 2 X 4 

+ XiX 3 )u 3 sin (2 gv — 2 0 O ) 

+ 6fc 3 (X 2 X 4 + X x X 3 ) 

•(1 — 5X x 2 X 4 2 )w 4 cos {gv — 0 O ) 

+ 6/c 3 X 1 X 4 (1 - 5X x 2 X 4 V 

•sin {gv — 0 O ) 

+ 30/b 3 X 1 2 X 4 2 (X 2 X 4 + X,X 3 )w 4 

• cos (3 gv — 30 o ) 

+ XOfeXiVu 4 sin {3gv — 30 o ) 

+ k,{3 - 60X x V + 210X 1 4 X 4 V 
+ /c 4 X x 2 X 4 2 (60 - 28 0X 1 2 X 4 2 )w 5 

• cos (2 gv — 2 0 O ) 

- fc 4 X 1 X 4 (X 2 X 4 + XjX 3 )(120 

- BCOXxW sin (2 gv - 20 o ) 

+ 70ft 4 Xi 4 X 4 V i cos (4 gv — 40 o ) 

- 280/c 4 Xj 3 X 4 3 (X 2 X 4 

+ XiX a )u 5 sin (4 gv — 4 0 O ) (17) 
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This development of the disturbing function 
is capable of producing the terms up to the 
second order in the mean motions of the argu- 
ments and the periodic terms of the first order. 
It can easily be extended if necessary. 
Perturbations in the orbit plane . We have 

u = h 2 + eh 2 cos (v — x) (18) 

h 2 = l/[o(l - e 2 )] (19) 

The orbital true longitude x of the osculating 
perigee reckoned from the departure point can 
be represented in the form 

X = (1 “ C)V + 7T 0 + <t> (20) 

where is a purely trigonometrical part which 
is of the order of perturbations; we obtain 

u = h 2 + h 2 e cos (cv — 7 r 0 — <£) (21) 

We introduce also 

u = h 2 + ho 2 e 0 cos (cv — 7 r 0 ) (22) 

where 

ho 2 = l/[do(l - eo 2 )] 

The equation 

l/u 2 dv/dz = l/K (23) 


But W this time is defined by the equation 
IF = — ^1 + ”^[1 + e 0 cos (cv — 7r 0 )] 

+ 2 — [1 + e cos (cv — 7 r 0 — <f>)] (31) 


or by the equation 



+ 2 ~ [1 + e cos (cv — 7 r 0 — <£)] (32) 

fi 0 

This is the fundamental form of Hansen’s W 
function if the orbital true longitude is used as 
the basic variable. 

The corresponding W function becomes 
W — — ^1 + “^(1 + e 0 cos <p) 

+ 2 — [1 + e cos (<p — <f>)\ (32') 

fi[) 

and 

w = w l.„.„ (32”) 

We have [Brown, 1896] 


which defines the pseudo-time z is equivalent to 
Kepler’s equation of the form 


and to 


Taking 


E — e 0 sin E = g 0 + cuqZ 


l/u = a 0 (l — e 0 cos E) 


l/u dv/dt = l/h 
into account we deduce from (23) 
dz/dt = h 0 /h u/u 


dbz 

dt 


h 0 u 


h u 


or the last equation can be put in the form 


1 * ho TF | 

j j - T 


dbz 

dt 


ho( v y 
h\l+J 


^ {he cos (0 — <p)l = — sin (/ + x — 0) 


(24) 


(25) 


+ 


dQ 

dv 


(r + p) 


cos (/ + x — 0) 


d& 


+ hef t sin ( X - 0), 


(33) 


(26) where j3 = (1 — c)v + 7 ro + <p. 
Taking 

1 + X = v 


(27) 


and, for the perturbations of time 8z = z — t, 
we obtain 


and (18) and (33) into account we obtain 

= h 0 -fr* cos ( cv — *0 — <p) 


dW 

dt 


(28) 


(29) 


— 1 — Co COS (p 

. 0 h 2 , , 

+ 2 7-5 { cos (ct> 
/lo 

2 dfl 


+ 


similar to the result in Hansen’s theory, where 
1 + v = u/u = r/f (30) 


h 0 dr 


To - <p) ~ ljj 

sin (cv — 7r 0 — <p) 


~ 2 ~ c ^ e 7t sin ^ - ^ ( 34 ) 
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Differentiating (32') with respect to <p, we have 

n h . ( v dW ( . h 0 \ 

— 2— e sin (<£ — </>)= + — Je Q sm <p 

and (34) takes the form 

, dQ f u , , 

= ho — 2 7-2 cos (cv — 7T 0 — ip) 

dv L 


dW 

dt 


— 1 — e 0 cos v? 


h 2 

+ 2 — 2 { COS (CV — 7T 0 — v?) — 1 ^ 
h o - 


, 2 dtt . , , 

+ ~r -r- sin {cv — ir 0 — </>) 

ho OT 


+ 

Taking 


.^-( 1 + ^) ei ) Sin ^ (1 " c ) l (35) 


dt _ A_L 

"T" “ ^0 * 1 2 

av n 0 u 


(350 


into consideration, we deduce the final form for 
the derivative of W to be used in this theory 


dW 

dv 


h 2 (- — — \ 9 — 

h ° w dv ) L 2 /io 2 


1 — e 0 cos <p 


COS ( CV — 7T 0 — (p) 


+ (1 - c) \% - 1 1 + f 2 )* 0 sb *] (36) 


+ 2 7*2 { cos ( cv — 7r 0 
too 

n d$l h . , 

— 2 — — sin {cv — 7To 
ow ho 


This form of the differential equation for W 
is much simpler than in the theory based on the 
use of the eccentric anomaly, and the divisor 
Vl - e„‘ does not appear. 

The integration of (36) leads to the equation 
of the form 


T = 2 — e cos ^ 1 + (39) 

= 2 — e sin <j> (40) 

ho 

Equations 38 and 39 contain the constants of 
integration in the left-hand side. These constants 
are determined in such a way that neither con- 
stant nor secular terms are present in no5z. From 
(38) we deduce, similarly as before, 



For the disturbed radius vector we deduce from 
(31) a new form for u: 

u = * + l{h n -~ 1 ) u + l h ° w i 0 (43) 


or 



The system of equations 36 and 41-43 is 
basic in the solving of the problem either by the 
numerical process of iteration or analytically. 
The numerical process is easier to produce the 
perturbations of any desired order. At each 
step of the iteration we will have the perturba- 
tions of different orders in the form of finite 
polynomials, as the development (17) shows, and 
not in the form of truncated trigonometrical 
series. 

If the analytical approach is taken, the results 
are obtained in the form 

+ k t 2 T t + k 2 3 T 3 + • • • 


W = S + T cos <p + 'k sin <p (37) where 


which is similar to equation 36 of the previous 
article, 


H = -I - 


t + 2 y 

h ho 


(38) 


Ti } T 2} 7*3, * * • 

are trigonometric polynomials in cv — tto and 
gv — 0o. Brenner [1960] was the first who recog- 
nized the existence of this development. 
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We deduce from (43) 

— = (l + — w) 7-^ (44) 

h \ u J 1 — v 

and taking (26) and (29) into consideration we 
obtain 



and 


dt hp 1 — v 2 

dv u 1 + hpW/u 


(46) 


it follows from the last two equations that 


drip hz 
dv 


_/■, „ 2 \ 3/2 

= (1 ” «o ) 



hp W /u + v 2 
* 1 + hp 2 W/u 


(47) 


or, in the form more convenient for the use of 
the process of iteration and with accuracy 
compatible with our information about the 
form of the earth, 


dn 0 5z 
dv 


= (1 - e 0 2 ) 3/2 (— V W 

\po/ 

■ n „ 2x3/2/ f V (« - uf 


_z w dn^i (48) 

Po dv 

We see that in the case of a large eccentricity 
the slowly convergent series of the form [Hansen, 
1855] 

M- _ 0 + jr 
W (1 -0 T' 1 

•|t> 0 <n> + 52 2 Vi (n) cos i (cv — 7T 0 )| (49) 

where 


So 


1 + y/ 1 — eo* 




(n) 


(-‘)v(” + ; _1 ) 


only in^the very last equation, without inter- 
fering with all intermediate results. 

The term of the form 

A + B cos {cv — 7 r 0 ) 

in (48) will disappear after the proper determina- 
tion of the constants of integration in IF. Only 
the first few terms in the series for 


{f/VoY (n = 1, 3, 4) (50) 

are important for the separating of this term, 
because of the form of the integrals. The con- 
vergence of the remaining parts of the series for 
(50) can be speeded up by a summability process. 

Perturbations of the orbit plane . Putting, as 
before, 

<r = (1 - g)v + 0 O - N - K (51) 

6 = (1 - h')v + Op - N + K (52) 

\ x = sin \i cos N 
X 2 = sin \i sin N 
X 3 = cos \i sin K 
X 4 = cos \i cos K 

and taking into consideration that j + ( w ) of 
the previous article becomes gv — 6 0j we can 
transform the previous equations for Xi, X 2 , X !; X 4 
to the form 


dXi (hf_ 
dt + \ 


+ 9 


■>•5 + 5 


rh 


dQ 


2 7”* dt ' 2 " dZ 

[+X 4 cos {gv — Op) — X 3 sin {gv — 0 O )] 


dX, 

dt 


--(*^-‘>•5 + 1 * 


dU 

dZ 


- [ X 3 cos {gv — 8p) — X 4 sin {gv — 0 O )] 


dX 3 

dt 


dX 4 


h' 


ZL 1 \ . 1 rh ^ 

2 X ‘ dt + 2 h 3Z 


■ [+X 2 cos ( gv — do) + X, sin (gv — 0 O )] 


h' ~ 9 , *>.1^30 

dt 2 3 dt 2 dZ 

• [— X, cos (gv — 6 0 ) + X 2 sin (gv — 0 O )] 

Taking (16) and (16')-(16'"), also 


•F(—n + 1, — n + i + 1, 1 + i, fi 2 ) 
reappear in the solution of our problem, but 


dQ _ 1 dtt 
dZ~ rd$ 
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into account, we deduce 

dX. _ | (V + a _ X * 

at + \ 2 1 r 2 <it 


, h 0 h dU 

+ T ho Mo 0051 


d \ o _ ( hf_ + <7 V dy 

dt “ V 2 / Al dJ 


h Q h dQ 

- T'jT^T cos 1 
4 /to oXi 


dX 3 


dX 4 


= 4 


h' — g x dy 


X 4 


di 

. /fc 0 h 012 

+ T'lT^T cos 1 
4 /i-o 0X 4 


h f — g dv 
~dt " " 2 Aa dJ 

Xo ^ 612 
4 /to 0X3 

Multiplying (54) by (35) we obtain 


cos i 


dX 

dv 


? - +(^ - >>* 


, a 0 2 ii ! iaa 

H 7“*T"2 “2 77“ COS 1 

4 /to it 0X2 


d\2 

dv 


- -(^ - *>■ 


feo 2 a 2 1 an 

4. '/l 0 2 7 dX, 


d\ 3 , h' — g 

dv + 2 X * 


, a 0 2 a 2 1 as) 

+ T'Tl 2 TT“ COS t 

4 /t 0 M 0X4 


dX, 

dy 


/j / ~ g 


A 0 2 A 2 1 dQ 

A ' 1 2 2 C0S l ~ 

4 h 0 u 0X3 


application of the process of numerical iteration 
as well as for the analytical development. This 
system (55) is incomparably more elegant than 
the system (31) of the previous article. The 
quantities h* and g are determined in such a 
way that no constant terms appear on rightr 
hand sides of (55). The position vector of the 
satellite can be written in the form 


(54) 


T ~~ u ^[0] • ^.1 [i] * A 3 [ a] 


where 


+ cos y 
+ sin v 
0 



+ cos a 

— sin a 

0 

A 3 [a] = 

+ sin a 

+ cos a 

0 


L 0 


0 

+ 1_ 


>1 


0 

0 

A, [a] = 

0 

+ 

cos a — 

sin a 


_ 0 

+ 

sin a + 

cos a_ 


We put 

(0) = (1 - a> + e u 

then 

0 = (0) _ n + K 
We also have 


(56) 

(57) 

(58) 

(59) 

(60) 

(61) 


<r = (1 — g) + 0 O ~ N - K 
Taking (59)-(61) into account, we deduce 

(55) r = i A 3 [(0)]-A 3 [K - N]- AS] 

cos ( gv — 0 O ) 

I 

■A 3 [K + N] 


sin (gv — 0 O ) 
0 


or 


r = - A 3 [(0)]-A- 
u 


The form (55) of the differential equations for 

Xi, X 2> X 3 , X 4 is new, and it is convenient for the where 


cos ( gv — d 0 ) 
sin (gv — 0 O ) 
0 


A = 


+Xi X 2 X 3 + X 4 2 , — 2(X 3 X 4 -j- X x X 2 ), +2 (XiX 3 — X 2 X 4 ) 

+2(X 3 X 4 XjX 2 ) , Xi "4“ X 2 X 3 H - X 4 , — 2(XiX 4 -f* X 2 X 3 ) 

L+2(XiX 3 + X 2 X 4 ), +2(XiX 4 — X 2 X 3 ), —X, 2 — X 2 2 + X 3 2 + X 4 2 J 
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A a[(0)] is the matrix connected with the regression 
of the node, and A gives small oscillations of 
the orbit plane around the mean position of 
the node. 

Conclusion. The theory of artificial satellites 
in terms of the orbital true longitude as given in 
this article has a simpler mathematical formula- 
tion than the theory based on the use of the 
eccentric anomaly. All formulas are more com- 
pact and have a wider field of validity. 

This theory is also more convenient from the 
standpoint of numerical computation. The 
programming will be simpler than for the 
previous theory, requiring less space in the 
memory, thus making the program more flexible 
and providing room for future refinements as 
our knowldege about the external gravitational 
field of the earth becomes more accurate. 
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